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ABSTRACT 

Late-time near-infrared and optical spectra of the normal-bright Type la supernova 2003du about 
300 days after the explosion are presented. At this late epoch, the emission profiles of well isolated 
[Fe II] lines (in particular that of the strong 1.644 /im feature) trace out the global kinematic dis- 
tribution of radioactive material in the expanding supernova ejecta. In SN 2003du, the 1.644 ^m 
[Fe II] line seems to show a flat-topped, profile, indicative of a thick but hollow-centered expanding 
shell, rather than a strongly-peaked profile that would be expected from a "center-filled" distribution. 
Based on detailed models for exploding Chandrasekhar mass white dwarfs, we show that the feature 
is consistent with spherical explosion models. Our model predicts central region of non-radioactive 
electron capture elements up to 2500-3000 km s~^ as a consequence of burning under high density, 
and an extended region of radioactive ^^Ni up to 9,000-10,000 km s~^. Furthermore our analysis in- 
dicates that the 1.644 /im [Fe II] line profile is not consistent with strong mixing between the regions 
of electron-capture isotopes and the ^®Ni layers as is predicted by detailed 3D models for nuclear 
deflagration fronts. We discuss the possibility that the flat-topped proflle could be produced as a 
result of an infrared catastrophe and conclude that such an explanation is unlikely. We discuss the 
limitations of our analysis and place our results into context by comparison with constraints on the 
distribution of radioactive ^^Ni in other SNe la and briefly discuss the potential implications of our 
result for the use of SNe la as cosmological standard candles. 
Subject headings: supernovae - nuclear burning fronts - progenitor systems 



1. INTRODUCTION 

There is general agreement that Type la Supernovae 
(SNe la) are the result of a thermonuclear explosion of 
a degenerate C/0 white dwarf (WD) most likely close 
to the Chandrasekhar limit for the majority of events, 
because these allow to reproduce optical/infrared light 
curves (LC) and spectra of SNe la reasonably well. How- 
ever, despite the success of Chandrasekhar mass mod- 
els, questions remain concerning the detailed physics of 
the burning front, the ignition process, the role of pre- 
conditioning and progenitor of the WD prior to the ex - 
plosive phase . For recent reviews see JRra.Tichl l|T999j) . 
IHoflich et all ll200l . iNoinoto et all Emm . 

The distribution of the iron-group elements is key to 
answer these questions. It has long been recognized that 
late-time spectra of SNe la, dominated by forbidden 
lines of Fe-peak elements, are an excellent pro be of the 
amount and distribution of radioactive ^^Ni ijAxelrodI 
Il980albt iMeverottilTflSnl) . After about day 200, the 
energy input is dominated by local energy deposition 
from positrons and the envelope is almost transpar- 
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ent, allowing for the direct observation of the entire 
distribution of radioactive matter, and making the 
result rather insensitive to radiation transport effects. 
However, the total amount of Ni is still rather uncer- 
tain because the emission depends sensitively on the 
redistri bution of energy within an ion (see also IR catas- 
troph e; IAxelrodlll980at iFransson. L unauist &: ChevalieJ 
1996ft. and o n the ionization mechanism (e.g . 

Spvromilio. Mei kle. Allen. & GrahanJ 1199^ 

Liu. .lefferv. & Schulta.1997; .Bowers et al.lll997f) . 

In principle, line profiles are less sensitive to these ef- 
fects since they depend not on the absolute amount of 
excitation at a certain time but rather on the relative 
change of the conditions within the envelope. Unfortu- 
nately the spectral features in the optical are produced 
by a large number of overlapping lines, making the in- 
trinsic emission proflles difficult to untangle. The near- 
IR contains much more isolated lines, potentially offer- 
ing the opportunity to look for the remnants of high- 
velocity plumes in the outer regions of the ^^Ni-rich gas 
and to look for kinematic offsets of the ^^Ni distribution 
from off-center ignition. However, infrared spectroscopic 
observations at such late epoch are difficult to obtain 
and to-date, very late-time (> 200d) IR spectra have 
only been publish ed for the unusual supernova SN 1991T 
Phillios ( im m-. iSnvromilio. Meikle. Allen, fc GrahamI 
l|1992il : iBowers et all 

Here we present a study of late-time spectra of the nor- 
mal bright Type la SN 2003du, concentrating on NIR 
spectroscopic observations with the 8.2m Subaru tele- 
scope. In § 2, we discuss the observations and data re- 
duction. In § 3, we describe our models used for the 
analysis of observations and its implications which is pre- 
sented in § 4. In § 5, we discuss the results with an eye 
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toward the broader context. In § 6, we sum up the re- 
sults and and comment on the hmits of our study and 
future perspectives. 

2. OBSERVATIONS AND DATA REDUCTION 

SN 2003du was discovered in U GC 9391 by LO- 
TOSS f Schw artz fc HolvorcemI l2?)?)3'l on 22 April 2003 
(UT) at a bout 15.9™. It was classified as Type la by 
iKotak fc M ciklc ( 2003) on 24 April 2003 who noticed 
that it resembled SN 2002bo about two weeks before 
maximum light. Early-epoch optical spectroscopy re- 
vealed a high-velocity component of the Ca II infrared 
triplet feature, similar to that seen in SN 2001el, which 
may be the signature of interaction with circumstellar 
material , but otherwise the optical evolution appeared 
normal l)Gerardv et alJl2004 - 

2.1. NIR Spectroscopy 

A near-infrared spectrum (1.1-1.8 /im R« 400, Fig.^ 
of SN 2003du, was obtained on 27 Feb 2004 (296 days 
after Vm.„,r) using the OH-airglow Suppressor (OHS 
llwamuro et al.ll2001|) and the Cooled Infrared S pectro- 
graph and Camera for OHS (CISCO Motohara et al.l 
12002) on the Subaru 8.2 m telescope. The observations 
consisted of two 2000 s exposures, using a 0'.'5 slit, with 
the object dithered 10" along the slit between exposures. 
A second, lower-S/N spectrum, totaling 1500 s on-target 
integration time, was obtained on 2 Apr 2004, (331 days 
after Vmax)- 

Data reduction was performed using standard IRAF 
routines fTody 1986). The two exposures were dif- 
ferenced for first-order sky subtraction, then the "neg- 
ative" spectrum was inverted, shifted, and combined 
with the "positive" spectrum. Cosmic rays, hot pix- 
els and other detecto r defec ts were remo ved using the 
LA-Cosmic package ijvan Dokkum 1200111 . Since the 
0'.'105 pixels oversample the spectrum, the combined data 
were boxcar filtered 3x3 before 1-D spectra were ex- 
tracted. During the extraction, a second background 
subtraction was performed using the average background 
level from the regions adjacent to the extraction aper- 
ture. 1-D spectra were wavelength calibrated using 
the standard OHS/CISCO wavelength mapping*, shifted 
to the measured slit position during the observations. 
The combined atmospheric and instrumental transmis- 
sion function was removed from the data by following 
the target observations with observations of an A2 star 
(SAO029603) at similar airmass, and matching the ob- 
served spectrum of the telluric star with stellar atmo- 
sphere models of Kurucz (1994). 

The S /N of the resulting spectrum is highly wavelength 
dependent. The data are completely sky- noise domi- 
nated so the S/N is essentially independent of the signal 
strength from the supernova, but is instead strongly de- 
pendent on the level of the night-sky emission at a given 
wavelength. Where strong night-sky emission remains af- 
ter passing through the OH-suppressor, a narrow region 
of low-S /N is created. This effect can be seen very clearly 
in the 2-D pre-extraction image shown in Figure^ where 
the bright O2 band near 1.3 /im is not suppressed by the 
OHS and creates a band of very low S/N in the resulting 
data. 

^ For details, refer to the Subaru OHS/CISCO websit e: 
|http: / / www.naoj .org/Observing/Instruments/OHS/index.html| 



The reduced NIR spectrum of SN 2003du on 27 
Feb is shown in Figure ^ with the wavelengths pre- 
sented in the rest-fra me of the host galaxy (cz = 
1914 km s-i; Schn eider erall Il992() . The spec- 
trum exhibits broad emission from the [Fe II] 
1.257 A*ni (a6D9/2-a'*D7/2) and 1.644 /^m (a^Fg/j- 
a^D7/2) lines which typically dominate nebular Fe 
spectra in the NIR (e. g. 'Nussbaumer fc StorevI '1980'; 
lOliva. Moorwood. fc Dan zigcr 1989. lj,9(i). Surprisingly, 
the 1.644 jjLm line appears to show a flat-topped line pro- 
file suggestive of an extended envelope with a central hole 
in the emissivity. In contrast, a filled envelope would 
produce a strongly peaked line profile. The S/N of the 
1.257 /im feature is too low to make any strong state- 
ments about the shape of the line profile, but it appears 
to be consistent with the boxy shape of the 1.644 /im 
feature. Similarly, the low S/N spectrum on 2 Apr is 
consistent with no significant evolution from the 27 Feb 
spectrum. Understanding the nature and implications of 
this apparent flat-topped emission profile is the subject 
of our analysis. 

2.2. Optical Spectroscopy 

In addition to the NIR observations, optical spectra 
(4200-8000 A, Rw 450) of SN 2003du were obtained on 
25 Mar 2004 (323 days after Vmax)i using the Imag- 
ing Grism Instrument on the 2.7 m Harlan J. Smith 
Telescope at McDonald Observatory. Observations of 
SN 2003du consisted of four 3600 s exposures using 
a 2" E-W oriented slit. Wavelength calibration was 
achieved by observing Cd, Ar, and Ne arc lamps and 
the sp ectra were flux calibrated using the Masscv et a| 
fl988) standards G191B2B, Hiltner 600, Feige 34, and 
PC 0823-^546. 

The reduced optical spectrum is shown in Figure 13 
The spectrum is typical of lat e-time optical spectra of 
SNe la (e.g. Bowers ct aLl ll997D . showing a strong blend 
of forbidden lines from Fe-peak elements in the 4000- 
5000 A region. The feature near 7300 A is the compar- 
atively isolated [Fe II] 7155 A fine, blended with [Ca II] 
7291 & 7324 A emission and other weak [Fe II] lines at 
7142, 7438 & 7648 A. 

3. MODEL CONSTRUCTION 

Our quantitative comparison is based on the delayed- 
detonation model 5p0z22.23 which reproduces the light- 
curve and spect ra of a normal-bright SN la like SN20 03du 
reasonable well l|Hoflich et alJ2002t[Ge'rardv et alj"2 004). 
In delayed-detonation models (DD), the white dwarf 
progenitor undergoes supersonic detonation after a pe- 
riod of pre-expansion , perhaps due to a phase of 
subsonic deflagration llKhokhlovl Il991t lYamaoka et alJ 
Il992t iWooslev fc Weaved 1199*^ We choose a de- 
layed detonation model because this class of model 
have been found to reproduce the majority of the ob- 
served optical/infrared light curves (LC) and spectral 
features of SNe la inc luding the brightn e ss decline re- 
lation (Hoflich 19 91 iFisher et alJ Il995t iHoflich et alJ 
19961 I' Hoflich & Khokhlov '1996^ ' Wheeler et all 11998 
Lentz et al. 2001; Hoflich ct al. 2002). As a reference, we 
also examine the profile i n the context of the well-known 
deflag ration model W7 l)Nomoto. Thielemann, fc Yokoil 
IT9M1 . 
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Fig. 1. — Near IR spectrum of SN2003du on February 27th, 2004, obtained at the Subaru telescope with the OH-Suppressor spectrograph 
(OHS). The 2-D data are shown after the data have been pair-subtracted and stacked, but prior to extraction and correction for instrumental 
response. Note that the S/N is highly wavelength dependent due to the influence of night-sky emission features even with the OH-suppressor. 
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Fig. 2.— Optical spectrum of SN2003du on March 27, 2004, normalized at 7000 A. 
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Our results are more generally applicable to all 
Chandrasekhar-mass models since the structure of the 
WD, the explosion energy, and the light-curves are 
mainly determined by nuclear physics rather than 
the details of the nuclear burning ( "stellar amnesia" ; 
IHoflich et alllaTnl . We use these specific models as 
testbeds to interpret the late-time observations. 

3.1. Open Questions of the Explosion Physics 

Despite the successes of delayed detona- 
tion models (and W7), questions remain con- 
cerning the detailed physic s of the b urning 
front ("Livne & Arnettl TgQl iKhokhloyI IToM 
2£inccke. Hillcbrandt. & Nieme ve^ Il999t IKhokhloyI 
2001lt . Rcincckc. Hillcbrandt . fc Niemeveil 2002; 
Gamezo. Khokhlov. ife Qran 20Q 4). the ignition process, 
and the role of p re-conditioning of the WD prior to the 
explosive phase ("Hoflich. Wheeler, fc Th ielemannl 119981: 
ITiscwski ct al. 2000: Hoflich fc Stein 200^. Overall, the 
propagation of the detonation front is well understood 
but the description of the deflagration front and, in 
the DD scenario, the mechanism for a deflagration to 
detonation transition (DDT) are persistent unsolved 
problems as is the location of the DDT which if off-center 
may result in off-center distributions of ^^Ni. 

On a microscopic scale, a deflagration propagates due 
to heat conduction by electrons. Though the laminar 
flame speed in SNe la is well known, the front has been 
found to be Rayleigh- Taylor (RT) u nstable increas ing the 
effective speed of the burning front ( Nom otclll982D . Cur- 
rent hydrodynamic calculations, starting from central ig- 
nitions in static WDs, all show that RT instabilities gov- 
ern the morphology of the burning front in the regime 
of linear instabilities, i.e. as long as perturbations re- 
main small. During the first seconds after the runaway, 
the increase of the flame surface due to RT instabilities 
remains small and the effective burning speed is rather 
small resulting in the rise of large plumes of burned mat- 
t er through the WD. 

iNiemever fc WooslevI ()1997|) studied the effect of off- 
center ignition and demonstrated that multiple-spot ig- 
nitions will alter the early propagation of the flame. In 
significantly off-center ignitions, the buildup time of RT- 
instabilities is shorter corresponding to the larger grav- 
itational acceleration, producing off-center distributions 
of iron-group elements, although the overall picture of 
large, rising plumes still remains unchanged. As a conse- 
quence (and in contrast to spherical explosion models), 
we expect mixing of electron-capture products, formed 
by high-density burning, with the radioactive ^^Ni-rich 
layers from lower-density burning. Note, however, that 
in all these calculations the constraint of starting from a 
static WD is not trivial, as the morphology (and, conse- 
quently, the effective burning speed) of blobs depends on 
small scale motions of the background produced prior to 
the explosive burning phase (Hoflich fc Stein 2002) par- 
ticularly in the earliest burning stage. 

Probing the properties of the deflagration burning 
front directly with observations has proven difficult be- 
cause the signatures of the deflagration are mostly wiped 
out during the detonation phase in normal-bright SNe la 
(jLivne 1999: Gamezo et al.l2 0'(13'l . This provides the jus- 
tification for our approach and the explanation why we 
can use a spherical DD model as baseline. 



3.2. Model Construction 

For this study, the delayed detonation model 5p0z22.23 
(Fig. O was chosen. We use also the deflagration model 
W7 as additional reference. 

We consider the explosion of a Chandrasekhar mass 
white dwarf which originates from a star with a main se- 
quence mass of 5 Mq and solar composition. At the time 
of the explosion, the central density is 2 x 10^ g cm~^. 
The nuclear burning starts as a deflagration with a pa- 
rameterize d descripti o n of r ate of burning based on 3-D 
models bv iKhokhlovl (I2001D . When the density reaches 
Ptr = 2.5 X 10^ g cm~'^, the detonation is triggered. The 
final density and chemical structures are given in Fig- 
ure 13 

The pre-expansion of the WD is determined only by 
the amount of burning during the deflagration phase and 
not the detailed properties of the deflagration front. In 
1-D spherical models ptr is thus a convenient parameter- 
ization of the amount of burning during the deflagration. 
After the detonation is triggered, the burning products 
depend mainly on the density under which burning oc- 
curs. Because the detonation wave propagates rapidly 
through the WD at slightly above the sound speed, the 
outcome is determined primarily by the density struc- 
ture of the (puffed up) WD at the time of the DDT. As 
a result, the flnal density and chemical profile is rather 
generic for all DD models. 

Our analysis of the line profiles is based on radia- 
tion transport models using our Monte-Carlo scheme im- 
plemented in HYDRA including detailed 7-ray trans- 
port and energy deposition of fast electrons (see 
Hoflich, 2003ab, and references therein). Forbidden 
lines arc included based on Kurucz's data base (1993) 
and Nussbaumer & Stor ey J 1988a b), and supple mented 
with line list s given in iLiu. .Tefferv. fc Schiilt j (^.997); 
IBowers et al.l lll997D. Guided bv detailed calcula- 
tions (jRuiz-Lapuente et al]ll995l: ILiu. Jefferv. fc Schultd 
|1997|) . we assume that the energy input is given by the 
energy deposition in form of 7-rays and positrons based 
on our detailed Monte Carlo calculations, and that the 
distribution of the dominant ion is given by the abun- 
dance profile in the envelope which allows us to inves- 
tigate line profiles. We do not calculate the ionization 
structure because its details depend sensitively on the as- 
sumptions and the processes. The total amount of Fe II 
varies and, as a consequence, the absolute fluxes are un- 
certain. For the Fe II distribution, the main difference 
between models in literature are related to the central 
region which may or may not show Fe II as the domi- 
nant ionization stage. Here, we use the observations to 
constrain the conditions at the center and, thus, have to 
use our approximation. 

Our approach is well supported by the current lit- 
erature. The most detailed calculations for the ion- 
ization structures near day 300 have been published 
for W7, a model with a chemical structure in the in- 
ner regions similar t o ou r delayed detonation m odel 
()Ruiz-LaDuente et alJll995HLiu. Jefferv. fc Schilit3ll997l 
hereafter RL95 and LJS97, respectively). Generally, the 
results of R95 and LJS97 roughly agree, namely Fe II is 
the dominant ionization stage throughout the ^^Ni region 
with little variation in the ionization fraction. Small dif- 
ferences, namely an increase of the Fe HI with radius, are 
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Fig. 3. — Chemical structure of the delayed detonation model. The abundances of stable isotopes and ^^Ni are given as a function of 
the expansion velocity. The curves with the highest abundance close to the center are ^''Fe, ^®Ni and ^^Fe. The calculation is based on a 
delayed detonation model with central and transition densities of 2 X lO'' g cm~^ and 2.5 X 10^ g cm"'', respectively. 



present at the outer, high velocity layers which contribute 
to the line wings. However, significant qualitative differ- 
ence can be seen close to the center where R95 find that 
iron group elements are mostly neutral and where LJS97 
find that Fe II dominates Fe I. At an epoch of around 300 
days, energy deposition and ionization by 7-rays is of lit- 
tle importance but collisions with fast electrons are the 
dominant source of ionization (LJS97). These fast elec- 
trons are produced by the energy deposition by positrons 
which are produced as a consequence of the f3~^ decay of 
^^Co. However, neither LJS97 nor R95 considered the 
effect of magnetic fields. Whether the fast leptons can 
leak into the core depends on the physical processes con- 
sidered and on the magnetic field or, more precisely, on 
the Gyro-radius R which is given by 

V2 me EeC 3 r 1 

R = ±2 3.4 X 10 cm 

eB± Bi_ ^ ^ 

with TOg, Ee and e being the mass, energy [in MeV] and 
charge of the fast electron, respectively, and B± the or- 
thogonal magnetic field. Typical progenitor WD mag- 
netic fields have be en estimated to be ~ 10^ - 10^ Gauss 
l)Leibundgutlll995l) . Assuming that the flux is frozen in, 
the magnetic field evolves as Br = Br{t = 0) x [ro/rY- 
Thus by day 300, the magnetic field is about ^ 10"^^ 
- 10~^ gauss, implying that the Gyro radius is about 
2% or less of the distance to the inner edge of the ^®Ni 
distribution (see Fig. |SJ). Even a small B± would keep 
the electrons local. For comparing the line profiles as a 
baseline model, we assumed that the local energy input 
is given by the ^^Ni distribution and that the emissivity 
is given by the abundance distribution. 

4. DATA ANALYSIS 

4.1. Delayed detonation scenario 

We will to concentrate on the spectral emission pro- 
file of the [FeH] at 1.644 /im which is essentially un- 
blended at late times. This feature is surrounded 
by weaker lines from the same multiplet, but in neb- 
ular conditions, these lin es emit at ~l/3 or less of 
the 1.644 /im feature (e.g. iNussbaumer fc StorevI 119801 



lOliva. Moorwood. & Danzigeill989tll990D . As such they 
may contribute some to the faint wings, but will not af- 
fect the bulk of the line profile. 

The comparison between the observation on 27 Feb 
and the models is given in Fig. ^ The flat-topped profile 
can be well reproduced by the theoretical line profile of 
our (unmixed) reference model (solid/blue line). In this 
model, the flat top is produced by the lack of radioactive 
matter (i.e. ^^Ni/Co) close the center as a consequence 
of high-density burning up to an expansion velocity of 
~ 3000 km s~^, resulting in non-radioactive electron- 
capture products (see Fig. O and also Brachwitz ct aQ 
l2000f) . The extended wings of the feature are produced 
by the ^^Ni/Co region which extends from ^ 3000 km s~^ 
to ^ 9, 000-10,000 km s^^. The slight asymmetry in the 
theoretical spectrum is caused by optical depth effects 
due to Thomson, bound-free and free-free opacities and 
is consistent with the observed profile, although the S/N 
is not nearly sufhcient to actually detect such a feature. 
In comparing with the observations, we shifted the model 
line profile 500 km s~^ to the blue. Such a shift might 
potentially be a signature of an off-center '^^Ni distribu- 
tion, but the shift is too small (w 1 spectral resolution 
element) to make any significant claim, especially given 
the S/N of the observed line profile. Due to the shorter 
integration time of the spectrum obtained on April 2nd, 
the level of noise is somewhat higher (Fig. [SJ. The IR 
feature at 1.644 /xm is still present in this spectrum at 
roughly the same flux level, and the profile is consistent 
with that seen in the earlier spectrum. 

In light of Rayleigh- Taylor instabilities predicted by 3- 
D deflagration models, the central concentration of elec- 
tron capture elements in our 1-D model might be re- 
garded as somewhat artificial. 3-D deflagration models 
produce a mixture of radioactive and non-radioactive iso- 
topes within the layers. To simulate this, we also calcu- 
lated the results from a model with mixing of the inner 
layers below 8000 km s~^. Similarly, a mixed model will 
also simulate the effect of penetration of fast electrons 
into the core in absence of magnetic fields. 

In Figure 01 the resulting line profile is shown in red. 
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Fig. 4. — NIR spectrum of SN 2003du on Febr. 27, 2004, in comparison with the theoretical line profiles (solid thick line). In addition, 
the individual components of the forbidden [Fe II] transition at 1.644 fini are given for the original delayed detonation model (solid) and 
mixed chemistry (light) normalized to the maximum line flux (dotted) and the wings (dashed), respectively. In addition, weaker individual 
contributions are given for weaker [Fe II] lines (1.599, 1.664 & 1.677 A). The theoretical spectrum has been shifted by the peculiar velocity 
of the host-galaxy (-|-1932km s~^) and a residual velocity of —500km s"'^. 



We show the profile with two normaUzations, one nor- 
malized to match the observed maximum flux and one 
adjusted to match the line wings. Although the S/N of 
the features is low, the observations do not appear to 
be a good match to either normalization of the mixed 
profile. The profile matched to the central emission is 
too narrow, consistently under-predicting the flux in the 
outer parts of the profile, while the fit to the wings pre- 
dicts significantly more flux in the center of the line pro- 
file than is observed. Since the S/N of the spectrum 
is admittedly low, it is difficult to absolutely rule out 
a more centrally peaked line profile, but the unmixed 
model certainly seems to be a much better fit to the 
data. Although the low-S/N line profile does not ex- 
clude all structures from rising RT plumes, the apparent 
boxy profile suggests that the rise of such plumes may 
be significantly limited. 

Given implications of the [Fe II] line profile, it is worth- 
while to test to see if the profile is present in other spec- 
tral features. Unfortunately, the S/N in the 1.257 /im 
feature is not sufficient to place much of a constraint on 
the line profile, though it is not inconsistent with the 
1.644 /im feature. In the optical spectrum, the features 



appear more peaked, seemingly in contradiction to the 
IR feature. However, the features in the 4000-5000 A 
are composed of blends of large numbers of forbidden 
lines from iron-group elements. Thus the profiles of these 
features are rather weakly linked to the kinematic distri- 
bution of radioactive material and, in our model at day 
300, do not allow to analyze its boxiness ^. 

The only relatively clean [Fe II] feature (7155 A) seems 
to contradict the observed 1.644 fim profile, showing 
what seems to be a definite peak at 7200 A(Fig|Hl). In- 
deed such a profile seems more suggestive of our mixed 
model (red) than the reference model with the central 
hole in the ^^Ni distribution (blue). However, neither 
interpretation explains the extended red wing of the ob- 
served emission. Instead, both the peaked emission at 
7200 A and the extended red wing extending to about 
7550 A can be well reproduced by two blended flat- 
topped profiles, consistent with the observed NIR line, 
produced by [Fe II] 7155 A and blends of weaker [Fe II] 
7438 A and [Ca II] 7291 & 7324 A which have an intrinsic 

Note that the profile will change with time and we did not 
study the time-evolution or a wide range of models 
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Fig. 5. — NIR spectrum of SN 2003du on April 2nd, 2004, with 1/3 of the integration time of the 27 Feb. spectrum. The hue is present 
and the profile is consistent with little evolution, still seeming to show a flattened profile. However, the S/N of the spectrum is such that 
it cannot put strong constraints on the shape of the line profile. The comparison line profiles are the same as for Fig. Bl 



emission ratio of 3:2. Note that the relative contribution 
to the red feature depends on the Ca abundance of the 
specific model. Some support in favor of a significant 
Call contribution may be taken from the observations 
because a sufficiently strong Fe feature at 7300 A would 
produce a secondary, red feature which has not been ob- 
served in SN 2003du (see Fig. EJ. The combined blend 
explains the red wing of the 'plateau' as the edge of the 
flat-topped [Ca II] line and, at the same time produced 
the steep red wing and flatter blue wing of the "[Fe II]" 
line through blends. Nevertheless, such ambiguity in the 
optical demonstrates the need of Near-IR observations at 
late time. 

4.2. Deflagration scenario 

Based on the spherical deflagration model W7, a sim- 
ilar analysis of the data has been carried out by Maeda 
using the code of Ruiz-Lapuente & Lucy (1992). In 
Fig. [7| we show a comparison of the observed [Fe II] 
1.644 /Ltm feature and a calculated line profile based 
on the deflagration m odels W 7 with some variations 
iNomot o. Thielemann, fc Yokoil fl984). We show the re- 
sults for the original, unmixed model (dotted) , with mix- 
ing the layers below 8000km s^^ (dashed), and assuming 
a hole in the emission for layers with velocities less than 
3000 km s~^ (solid) with the flux normalized to the peak 
of the unmixed model without adjusting the profiles for 
the other model to the observed data. At day 300 and 
the resulting profile is not entirely flat topped due to two 
effects: first in W7, ^^Ni can be found down to about 



2200 to 2300 km s^^ and, second, this calculation still 
exhibits significant energy deposition by 7-rays near the 
center. In part, the central 7-ray deposition is caused by 
the assumption that each scattered results in full ther- 
malization. Apparently, we require a slightly larger hole 
of 3000 km s~^ to create a sufficiently broad, flat topped 
feature. 

5. ALTERNATE INTERPRETATIONS 

Our suggestion that the flat-topped profiles of 
SN 2003du indicate a concentration of electron-capture 
isotopes in the central region with no or little mixing 
in of ^^Ni, is not the only possible interpretation. For 
examp le, in the 'infrared catastrophe' picture ijAxelrodl 
Il980ajl . a lack of emission at certain expansion velocities 
can be caused by a redistribution of energy from the opti- 
cal/near IR to the far IR when the density drops below a 
certain limit (about 100 days after the explosion). How- 
ever, the brightness of SN 2003du is consistent with the 
predictions from early times rather than showing a rapid 
drop. Furthermore, the outer, low-density, high-velocity 
layers should exhibit the effect first, and we see no evo- 
lution of the line profiles between Feb 27 and Apr 2nd, 
2004. Moreover, the optical spectrum looks typical for 
late-time spectra. Thus, as none of the other predicted 
signatures of an IR catastrophe can be identified in the 
observations, we regard this alternative explanation as 
unlikely. 

A difference in the ionization structure of iron from the 
one assumed here could be another explanation since the 
observed line profile depends on the Fe II distribution 
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Fig. 6. — Comparison of the observed and theoretical emission profiles (solid, dark grey) for the region between 6800 and 7790 A. Fitting 
the line fluxes, the black and Ught-grey lines (sohd) present fits of the [Fe II] (7155 & 7172A) and [Ca II] (7291 & 7324A), respectively. 
Contributions of [Fe II] lines are given for other weaker transitions of [Fc II] (dotted and dashed-dotted). The dotted line shows the 
[Fe II] (7155 & 7172A) for the mixed model. 
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Fig. 7. — A flat-topped profile can be produced by any model without emission for velocities < 3000 km s"'^. We show a 

comparison of the observed [Fe I I] 1.64 4 fim feature and a calculated line profile based on the defiagration model W7 with some variations 
INomoto. 'I'hielemann. fc Yokoij 119841) . We show the results for the original, unmixed model (dotted), with mixing the layers below 
8000km s-l (dashed), and assuming a hole in the emission for layers with velocities less than 3000 km s ^ (solid). The fluxes are 
normalized to the peak of the unmixed model without adjusting the proflles for the other model to the observed data. 
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and, in principle, a shift in the ionization state to either 
Fe I or III may produce flat-topped Fe II profiles. In prin- 
ciple, a shift to Fe III may show up in the region around 
4700 A. However, in practice, such analysis is severely 
hampered by strong blending in this region which pro- 
duce very wide features equivalent to « 20, 000 km s^^ 
compared to 3,000fcTO considered here. Deconvolv- 
ing such a feature from this blend depends sensitively on 
details of the ionization and excitation structure. In our 
observations, we do not see evidence for a strong narrow 
emission component by [Fe III] in this region though, 
empirically, some contribution cannot be ruled out. 

Detailed calculations for the ionization structures at 
about day 300 have been published for W7, a model with 
a chemical structure similar to our delayed detonation 
model (RL95; LJS97). In these models the ionization 
balance shifts from Fe II to Fe III from the inner to the 
outer regions because the recombination rate and charge 
exchange reactions decrease quadratically with the den- 
sity in the region with radioactive isotopes. As a conse- 
quence, higher ionization stages are seen at larger radii, 
although the Fe II/Fe III ratio is rather smoothly decreas- 
ing with density, without sudden changes in the ioniza- 
tion balance. Thus, the lack of Fe II emission from the 
central layers cannot be attributed to a change in the 
balance towards Fe III. Furthermore, in the radioactive 
region, ionization by (local) positron deposition keeps the 
Fe II/Fe I ratio high. In presence of mixing, we cannot 
expect a 'hole' in the Fe II distribution which leads us 
back to our basic scenario of a central hole in the ^^Ni 
distribution. 

The necessity of a hole of 3000 km s'^ in the ^^Ni 
distribution is not unique to these specific models we 
have examined here. Any Mch model will require a sim- 
ilar hole to be consistent with the observation. The size 
of the central region with elements produced by electron 
capture may, in principle, provide a measure for both the 
central densities and possible mixing. Currently however, 
uncertainties in the treatment of the deflagration front 
and the electron capture rate do not allow us to discrim- 
inate between specific models ( Brac hwitz ct al...2000.') . 

6. DISCUSSION 

6.1. SN la Explosion Physics 

Although the observed profile can be reproduced by 
spherical models of massive WDs, our results pose a 
serious challenge to our understanding of the detailed 
physics of the explosion. All current generation 3-D mod- 
els predict deflagration fronts with properties dominated 
by slowly-growing large-scale Rayleigh- Taylor instabil- 
ities which allow the products of high-density nuclear 
burning to rise in plumes and mix with ^^Ni. However, 
our results suggest that the early phase of the explosion, 
which leads to the early expansion of the WD, is not 
dominated by RT instabilities (at least not in the case of 
SN 2003du), though there were similar suggestions from 
the analysis of SN 1999by a very subluminous SNe la. 

Perhaps the early phases of nuclear burning are dom- 
inated by pre-existing veloci ty fields produced pr ior to 
the thermonuclear runaway fflofiic h fc Steh]l200^ . Al- 
ternatively, this phase might be dominated by the rise of 
a single large bubble of burned material which penetrates 
the WD surface and triggers a detonation that propa- 
gates inward (Plewa et al. 2004). We may speculate 



that variations in the central density of the WD may be 
responsible for the appearance of profiles like that appar- 
ently seen in SN 2003du or, even more radically, that the 
pre-expansion of the WD prior to the detonation phase is 
not produced during the explosive phase of burning but 
rather during the pre-condoning phase followed directly 
by an explosive detonation phase. 

6.2. Comparison With Other SNe la 

It is useful to put our findings into context and exam- 
ine the signatures of mixing the very inner layers in other 
SNe la. As mentioned above, the peculiar SN 1991T is 
the only other SNe la for which late-time (> 200 d) NIR 
line profiles have been published. In contrast to our ob- 
servations of SN 2003du, the line profiles in SN 1991T 
were triangular-shaped, strongly indicating mixing of the 
inner layers. Late-time observations of the "normal- 
bright" SN 1998bu (Meikle 2004, Spyromiho et al. 2004) 
also seem to show a more centrally peaked profile in the 
NIR [Fe II] lines 

Note that in the scenario we suggest here, even in a 
SN with very little mixing, the flat-topped profile will 
be a transient phenomenon. At earlier times, the central 
regions of the ejecta will still be optically thick to 7-rays, 
and will therefore show significant low-velocity emission. 
As the inner layers become optically thin to 7-rays, the 
central peak would be expected to fade, leaving a boxy 
profile like the one discussed here, if the magnetic field in 
the ejecta is strong enough to prevent the free streaming 
of fast electrons from the radioactive ejecta into the cen- 
ter region. Eventually the expansion of the ejecta would 
dilute the magnetic field sufficiently to allow the pene- 
tration of fast electrons. The energy deposition would 
no longer be local, and the central peak would light up 
again. Thus, this scenario would predict an evolution 
where the [Fe II] features start centrally peaked, proceed 
through a phase of boxy emission, and then return to a 
peaked emission phase (unless the magnetic field in the 
WD progenitor is too weak to cause the local trapping 
of fast electrons, in which case the profile would always 
remain peaked.) Indeed, if such an evolution would be 
observed in the late-time spectra of a Type la SN, the 
epoch of the second transition to a peaked profile could, 
in principle, be used to infer the strength of the magnetic 
field in the central region of the WD progenitor. 

The very subluminous SN 1999by presented a different 
constraint on the kinematics of the ^^Ni-rich ejecta. In 
subluminous SNe la, little ^^Ni is produced during the 
detonation and the Si-rich region extends down to expan- 
sion velocities of « 5000 km s~^. As a consequence, the 
chemical signatures of the deflagration phase will sur- 
vive the subsequent detonation and has been expected 
to show up as iron-rich plumes embedded in a Si/S rich 
region. However, strong mixing can be excluded based 
on the spectral evo lution o f SN 1999b y at early times 
l|HoweU et aLllMl Hoflich etHMoO^ . 

6.3. Implications for SN la Cosmology 

Our results not only pose a challenge to our present 
understanding of the physics of the explosion but also 

However, e.g. the optical spectrum of SN 2003cg around 
400 days after maximum light shows a very boxy shape to the 
[Fe II] /[Ca II] emission feature near 7200 A (Elias de la Rosa et 
al., 2004) 



10 



Hoflich et al. 




20 40 60 

time [d] 




time [d] 



Fig. 8. — Visual light curves and (B — V) for the same delayed detonation model but with and without mixing of the inner layers (expansion 
velocities < 8000km s~^). Our LC calculations take into account detailed hydrodynamics, nucl ear netw orks, radiation transport for low 
and high energy photons, opacities and include solvers to calculate the atomic level populations lHoflichll20031) . 



raise questions of whether SN 2GG3du is an exception 
or the norm, and how this will affect our understand- 
ing of the diversity of SNe la. The implications for 
the use of supernovae in high precision cosmology may 
be sig nificant as it requires a 2 % accuracy in dis- 
tance (jOst riker fc Steinhardt"200r: 'Wcllcr J. fc Albrecht' 
I2OOI . In Fig. El light curves arc shown for two explosion 
models (Fig. with and without full mixing of layers 
with expansion velocitie s < 8 000 km s~^ (guided by 3-D 
calculations (Khokhlov 1200 ll) *). but otherwise identical. 
Such mixing will bring some ^^Ni to the central, high 
density region and reduces the ^^Ni abundance at other 
layers. Quantitatively, the results will depend details of 
the mixing, and the differences are of the same order as 
the numerical accuracy of our calculation. Thus, the LCs 
are shown here to demonstrate the qualitative effects and 
the approximate size of the effects. 

The difference can be can be understood in terms of 
energy deposition by the 7-rays and the diffusion time 



scales of low energy photons in the SN envelope. The lat- 
ter st rongly decreases with time t (roughly c>c t~'^ .HoflichI 
11995ft and, at maximum light, it is comparable to the 
expansion time scale. Up until about 2 to 3 days be- 
fore maximum light, the luminosity is mostly determined 
by the outer part of the ^^Ni layer and, consequently, 
the early light curves are virtually identical. With time, 
deeper layers start to contribute. Because mixing effec- 
tively reduces the ^^Ni which contributes to the lumi- 
nosity prior to about day 30, the LC of the mixed model 
declines faster and is redder. A competing, opposite ef- 
fect is related to the energy deposition by 7-rays or, more 
precisely, the escape probability for high energy photons. 
Increasingly, 7-rays produced in the outer layers of the 
^^Ni distribution can escape and, as a consequence, are 
lost for heating the envelope. In SNe la, the escape prob- 
ability for 7-rays increases strongly with time from about 
10 to 15% at maximum light to about 60% at about day 
50 (Hoflich et al. 1992). Mixing of ''^^Ni increases the 
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trapping of 7-rays because the inner layers are optically 
thick. As a consequence, the himinosity of the mixed 
model is higher and the temperature slightly higher. In 
conclusion, mixing processes will affect the brightness de- 
cline relation and the peak to tail ratio on the 0.1 mag 
level. Depending on the amount of mixing and its de- 
pendence on the progenitor, this may cause small but 
systematic evolutionary effects with redshift. 

7. SUMMARY 

We have presented late-time near-infrared and optical 
spectra of the 'normal-bright' Type la supernova 2003du. 
We show that late-time infrared spectra provide a probe 
of the distribution of radioactive ^^Ni in the ejecta. The 
strong, non-blended [Fe II] emission line at 1.644 /im is 
particularly well suited for such work, and in SN 2003du 
it appears to exhibit a box-like profile, whereas the opti- 
cal [Fe II] features (including the strong, single [Fe II] line 
near 7200 A) appear more peaked like those previously 
seen in other SNe la at comparable epochs. 

Based on detailed models for exploding 
Chandrasekhar-mass white dwarfs, we have presented 
a quantitative analysis of the observations, including 
optical depth effects (although these are rather small). 
Using our detonation models, we have shown that the 
line profile of the 1.644 /im line is consistent with spher- 
ical explosion models which produce a central region 
electron-capture elements up to « 2500-3000 km s~^ as 
a consequence of burning under high density, and an ex- 
tended region of radioactive ^^Ni up to w 9, 000 10,000 
km s^^. We have shown that the observed IR line profile 
is inconsistent with mixing of the region of neutron-rich 
isotopes with the ^^Ni-rich layers. The 1.644 /xm line 
is offset about 500 km relative to the host galaxy 
which might indicate a slight shift of the distribution of 
radioactive elements with respect to the envelope. 

The observed peaked profile of the optical [Fe II] fea- 
ture at 7200 Aand its extended wing can be well repro- 
duced by a blend of [Fe II] 7155 & 7172 A and weaker 
fines around 7388 & 7452 A , and [Call] 7291/7324 A, 
with line profiles consistent with the 1.644 /im feature. 
Such ambiguity in the interpretation of the optical lines 
demonstrates the usefulness of near-IR observations be- 
yond day 250-300, at phases when 7-ray deposition be- 
comes imimportant for the ionization balance of iron 
group elements. 

Although a hole in the ^^Ni seems to be a neces- 
sary condition to produce a flat-topped line profile, 



we also need to require sufficiently local trapping of 
fast electrons, as could be produced by weak mag- 
netic fields. Higher signal-to-noise observations of other 
normal-bright SNe la, coupled with detailed ionization 
models are certainly needed and will provide insight into 
the diversity (and causes thereof) in SNe la, and also into 
possible sources of evolutionary effects at higher redshift. 

While the fiattened [Fe II] 1.644 /<m profile is consis- 
tent with our 1-D models, it presents a challenge to our 
understanding of the physics in the earliest phase of the 
explosion, since all current 3-D models of the deflagration 
phase predict significant mixing in the innermost regions. 
This result for SN 2003du, while differing significantly 
from the late-time observations of the unusual SN 1991T, 
places constraints on the mixing of ^^Ni during the de- 
flagration which are similar to constraints placed on the 
mixing in the extremely sub-luminous SN 1999by. We 
now have examples of SNe la in which mixing may be 
suggested (SN 1991T, SN 1998bu) and in which it is ex- 
cluded, or at least very much limited (SN 1999by and 
SN 2003du). Understanding these phenomena will have 
significant implications both for the physics of SN la ex- 
plosions, and for their use as cosmological probes. 

Finally, we also have to address the limitations of this 
study which are subject to ongoing projects. Programs 
at SUBARU will increase the number of SNe la with late- 
time IR-observations and to study the spectral evolution 
with time. Detailed calculations are under way for the 
ionization structure which include magnetic fields with 
the goal to study their influence on the spectra and the 
spectral evolution. 
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